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Abstract. - Magnetic relaxation in a type-II superconductor is simulated for a range of tem- 
peratures r in a simple model of 2D Josephson junction array (JJA) with flnite screening. The 
high-T phase, that is characterised by a single time scale t^, crosses over to an intermediate 
phase at a lower temperature Tcr wherein a second time scale rp <^Ta emerges. The relaxation 
in the time window set by follows power law which is attributed to self-organization of the 
magnetic flux during relaxation. Consequently, for T < Tcr, a transition from super-critical 
(current density J > Jc) to sub-critical (J < Jc) state separated by an intermediate state with 
frozen dynamics is observed. Both and diverges at Tsc < Tcr, marking the transition into 
a state with true persistent current. 



Introduction. - In a hard type-II superconductor, flux-creep over energy barrier at a finite T 
leads to magnetic relaxation over long time scales. The flux creep, occurring due to thermally 
activated hopping of the flux lines, tend to reduce the local field gradient and hence the 
current density J. In analyzing relaxation measurements, the initial magnetic field distribution 
is assumed to be that of Bean's critical state |l| wherein J is replaced by the critical current 
density J^,. If the relaxation is close to Jc, the magnetisation decay M{t) is then theoretically 
known to be logarithmic, as is the case in most low- Tc superconductors ||^. In this case, the 
effective pinning potential U{J) Uo{J — Jc) is a good approximation. 

In high-Tc superconductors (HTSC), the large thermal energy available leads to rapid 
decay of M(t). Experimentally the relaxation is observed to be non-logarithmic over several 
decades in time in these materials This is interpreted as arising from a non- linear U{J). 
The vortex-glass theory Q and collective-creep theory which predicts a low temperature 
true superconducting state with finite Jc, expects a pinning potential of the form U{J) = 
Uo[{Jc/ J)^ — 1] which diverges in the limit J —^ 0. An important experimental observation in 
HTSC is an apparent universal value of the normalized relaxation rate S{T) = -jj-j^ around 
0.02-0.035 over a wide range of T For the C/(J) given above, such a small T-independent 
S requires fi > 2 which is beyond the range of existing theoretical models Q. Experiments 
have also indicated power law decay in HTSC which requires a logarithmically diverging 
U{J) = Uol'n{Jc/ J) 1^, ^. This form of U{J) cannot account for the experimentally observed 
plateau in S{T). 
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It becomes then important to identify the relaxation behaviour of the magnetisation in a 
type-II superconductor in presence of a uniform background of pinning centers. Towards this 
end, we simulate the time decay of remanent magnetisation M{t) at finite T in a simple model 
of 2D Josephson junction array (JJA). The magnetisation in this model is studied by including 
the screening currents through the inductance Lr.r/ between the cells. The underlying discrete 
lattice of junctions provides an energy barrier for the vortex motion within the array ||ic| and 
is the source of vortex pinning in JJA. The behaviour of JJA with screening is parameterized 
by A J = 27rLo/, ' where Ic is the critical current of the junction, Lq is the self-inductance of the 
cell, and $o is the quantum of flux. Detailed simulation at T = have shown that for Aj < 1, 
the magnetic response of this model is similar to a continuum hard type-II superconductor 
[ pi] p^ , p^ , In this paper, we show that the thermo-remanent relaxation of M{t) of JJA 
captures essential features of magnetic relaxation in HTSC. At a crossover temperature Tcr, a 
new time scale Tp emerges in an intermediate time window which along with the characteristic 
time scale Ta for long time behaviour governs the relaxation. The characteristic time scales Ta 
and Tp diverges as a power law at a temperature Tsc at which Af (r) — s- Mq ^ as t — > oo. 

The Model. - We consider a 2D array of superconducting islands forming a homogeneous 
square lattice oi N x N unit cells in the x—y plane. Tunneling of the macroscopic wavefunction 
^(r) = Tp exp[i(p(r)] across neighbouring islands lead to Josephson coupling between them. In 
presence of an applied magnetic flux $ext (per cell) along z direction, the junction behaviour 
is fully described by dynamics of the gauge-invariant phase difference 0r,(5 = 'p(r) — <y5(r + 5) — 
^ J^^^ A ■ dl between neighbouring islands. Here, (5 is a unit vector, and A is the vector 
potential corresponding to the total magnetic field. The inset of Fig.l shows a schematic array 
of size N = 5 along with variable (j)r^x and (t)r,y 

The magnetic response of JJA is set by the induced flux (R) due to screening currents 
which is modeled by considering the geometrical inductance matrix L of the array [T^ . This 
allows us to write the induced flux in a cell at R as <I>i„d(R) — J2r' L{'R,Il')I{'R') where 
/(R') is the cell current at R' and i(R, R') is the mutual inductance between the cells at R 
and R'. In order to ease the prohibitive computation cost involved when mutual inductance is 
considered jl^, Q, we consider the induced flux only due to self-inductance Lq ~ i(R, R) of 
the cell. This approximation is equivalent to the case of a long (ideally infinite) superconductor 
parallel to an applied field for which demagnetisation factor Nu = 0. The total flux at R then 
can be written as $r — $ext + ^o^R (the cell co-ordinates are used as subscripts). The cell 
current is only a convenient variable for introducing $i„d as the current through the junction 
is 1-rfi = Ir — Ir-5 where r is the junction common to cells at R and R — J. Since, /r is 
divergence-less, Kirchoff 's law is automatically satisfied at each node of the lattice. 

The dynamical variable (jj^ (subscript S is implicit) is related to the cell current /r through 
the flux-quantisation condition 



where the summation is taken around the cell in anti-clockwise direction. The equation of 
motion for (j)j.^s then follows from 




where the Hamiltonian Ti. of the system i; 
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(3) 
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where T = (f^)^-^ with R and Ic as the normal state resistance and critical current of the 
junction respectively. In the above equation, the first term is summed over all bonds and 
represent Josephson coupling energy for the junction, whereas the second term is the magnetic 
field energy due to screening currents and the summation is over all cells in the array. 

The above equations can be written compactly by introducing the matrix M for the lattice 
curl operation Eq.(l) then becomes M(/)~~27rf—^LoI„i where (f> and I„i are the column 
vectors formed by 0r and /r respectively, and / = ^ext/^o- Also, the current through the 
junction Ii, — M"^/™, where M"^ is transpose of the matrix M. For finite temperature simulation, 
we couple the heat bath through the noise current in the shunt resistor R across the junction. 
Langevin equation for the array then takes a simple form 

sin0 + ^{t), 

^Ira- (4) 

Here, the current Im is scaled by the critical current Ic of the junction. The dimensionless 
time T = t, and Aj has been defined earlier. The random term X{t) has zero mean and 
white noise correlation 

(Xr(r))=0, and (Xr(T)Xr' (r')) = 2T5(t - T')^r,r' , (5) 

where T is the temperature of the bath (in units of Ic^o/2nkB)- In this unit, the phase 
coherence between neighbouring islands is established at T = 1 and is the superconducting 
transition. The set of equations in (4) is solved self-consistently at each time step with 
free-end boundary condition. For simplicity, Ic is assumed to be independent of / and T 
(the / dependence docs not show any qualitative change in the results presented below). The 
magnetisation M is obtained as AI ~ {l/N"^) X]r(*5^r/^o) — /• 

The simulations were performed for N — 16(= 256 cells) with time step Ar = 0.05. The 
results presented below are for Aj = 0.1. Note that for Aj < 1, each cell accommodates 
more than a quantum of flux. The moderate size of the array allowed long time to be 
reached in the relaxation that is essential at low temperatures (the longest relaxation time 
run was of « 5 x lO^r w 10^ iterations). The field cooling is done by quenching from a 
high temperature (typically T = 2) to the temperature T under consideration, followed by 
annealing for 5x10'^— SxlO'^r before switching off the applied field at r = (thermoremanent 
magnetisation) . We have also carried out simulation for slow cooling which shows no qualitative 
difference from that obtained after a sudden quench, though a small but discernible dependence 
on the cooling rate is observed. We defer such effects and further details to a future paper, 
and present here results which brings out generic features. 

Results and discussions. - 

The AI{t) is shown in Fig.l over 7 decades of time after cooling to different temperatures 
in an applied field f — b (curves for intermediate T is not shown in order to avoid over 
crowding of the figure). The curves are scaled by M(t = 0). From the curves, two distinct 
temperatures Tcr ~ 0.26 and Tsc ~ 0.04 can be identified at which the dynamical behaviour 
changes remarkably. For T > Tcr, M{t) can be characterized by a single time scale Tq,. At 
T = Tcr, M{t) develops a kink in an intermediate time window which at lower temperatures 
evolve into a plateau. The now characterizes the long time behaviour of M{t). On the 
plateau, the magnetisation M(t) — Mq is T-independent. The temporary freezing of dynamics 
indicates emergence of a new time scale, represented by rp, which governs the relaxation for 
T << Ta- Both Ta and Tp increases rapidly at lower temperatures, and at Tsc, Mir) freezes 
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asymptotically to Mq for the time probed in simulation (note that for Tgc < T, M{t) becomes 
zero in the limit r — > oo). The temperature T^c thus marks the transition into a state with 
true persistent current. For T < Tgc, M{t) relaxes towards Mq. 

Also shown in the inset of Fig. 1 is the magnetic susceptibility x — ^'^1 f under field cooled 
(FC) and zero-field cooled (ZFC) conditions (for ZFC, / is applied at T = and x is calculated 
while increasing T). For T > Tcr, the xfc and xzfc are equal and M{T) is reversible. At 
Tcr, the difference | xfc — Xzfc |> 0, and the xzfc below it depends strongly upon the 
thermal history [|l5j. Similar behaviour in bulk superconductors allow us to identify Tcr as 
the irreversibility temperature. Though xfc does not show any change at Tsc, Xzfc is T 
independent below it. The T independent value of xzfc ~ 0.6 gives Mzfc ~ 3.0 which 
equals remanent magnetisation Mq obtained from the relaxation for T < Tsc and is also found 
to hold for other values of /. This equality between the Mq and Mzfc can be understood 
by invoking the relation Mrem{H) = Mfc{H) — Mzfc{H) {H is the applied field) which is 
experimentally known to be valid in bulk type-II superconductors [^6| . In case of strong flux 
pinning, Mpc ~ and Mrcm ~ —MzfC: which is the case observed here. Validity of the above 
relation is also an evidence for strong flux pinning in JJA arising due to discrete underlying 
lattice of junctions. Further, we analyze the relaxation occurring on different time scales. 

The long time decay of M(r), which we term as a-relaxation, fits Kohlrausch- Williams- Watt 
(KWW) law M(r) = exp[— (;^)"] over the temperature range T > Tsc probed in the simulation. 
The Ta (T) shows 6 orders of increase between Tsc and Tcr implying that the system falls out of 
equilibrium on cooling through it. Though Ta(T) appears to fit Arrhenius law TalT) ~ e"^/^ 
(see Fig. 2 inset (a)), notably, a distinct power law is observed at low temperatures as shown 
in Fig. 2 where is plotted against T — Tsc- Fit to Ta{T) = tq \ T — Tsc I"'' gives 
Tsc = 0.045,7 = 3.74 for / = 5, and Tsc = 0.05,7 = 4.735 for / = 2 which is also included 
in the figure. The value of Tsc obtained from the fit is not very different from that at which 
the relaxation is frozen asymptotically at Mq but it must be treated with considerable reserve. 
It is probable that the actual value may be lower than that obtained here as the time scale 
probed at low temperatures put severe constraint on obtaining Tq. In Fig. 2, the deviation 
from linearity occurs at T w 0.24 which is close to Tcr below which the relaxation shows a 
temporary frozen state. The T dependence of the exponent a is also shown in Fig. 2 (inset (b)). 
Due to small value of a at low temperatures, the M(t) can be fit to logr over 1-2 decades in 
time, and is the regime where thermally activated flux-creep theory can be applied |^. 

Appearance of a new time scale Tp for T < Tcr leads to two-step relaxation process : an 
initial part towards the plateau during tq <^t<Ct^, and a later part away from the plateau 
in the interval Tp <^t <^Tci which at later time develops into a-relaxation. Here, tq « 20t is 
of the same order which appears in the power law fit for Tq,. We refer the relaxation during 
To <C r <IC Tq, as the /3-relaxation regime due to apparent qualitative similarity with the 
/3- relaxation seen in supercooled liquids p^ . For Tsc < T, the relaxation during Tp <^Ta 
fits a power law M(r) — Mq ~ —c{T/TpfoveT 3 decades with exponent b independent of T. 
This is shown in Fig. 3 where M(t) (for f — 5) for different temperatures can be scaled on to a 
master curve when plotted against T/Ta{T). The thick line is the power law fit with h = 0.85 
and holds for T/Ta{T) ^ 1 as evident from the fit. The scaling exponent b allow us to obtain 
the T dependence of Tp{T) which is shown in Fig. 3 (inset). Tp{T) fits to a power law of the 
form Tp = tq \ T — Tsc \~'^ with ij} = 1.58 and Tsc = 0.033 (tq is of the same order as that 
obtained for Tq,). Since the /3-relaxation at later time evolves into a-relaxation, divergence in 
Tp and Tq must occur at the same temperature. We attribute the small difference seen here 
to uncertainties in obtaining Tq. Nevertheless, divergence of Tp and Tq at Tsc is significant as 
it implies absence of flux-creep (hence, divergent U{J)) as t —>■ go, and is consistent with the 
idea of true superconducting state below a transition temperature. 
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For to^t<Ct^, the relaxation occur towards the plateau, and for T < Tsc is asymptotically 
frozen on it. The relaxation can again be fit to a power law M(t) — Mq ~ c' {t / Tp)~°' . The 
exponent a is dependent on T as shown in Fig. 4. Also shown in the figure is S{T) close to 
plateau which attains a T independent value « 0.076 (on the plateau, S{T) is an order of 
magnitude lower than this). This value is of the same order as that observed in HTSC for 
which it falls in the range 0.02 — 0.035 A plateau is S{T) coupled with non-logarithmic 
decay of M{t) over 4 decades is significant in view of similar experimental observation in 
HTSC. To understand the processes involved during the relaxation, the spatial distribution 
of $R, in different time windows is obtained and is shown in the inset of Fig. 4. The slope 
of the flux profile (d^/dR^) cx J is observed to be T- independent on the plateau and equals 
the slope (oc Jc) in the remanent state at T = (note that Ic of the junction is assumed 
to be T-independent in the simulation). At T = 0, such a state have been shown to be the 
self-organized Bean's critical state for a hard superconductor |l2). 

Thus, across the /3-relaxation regime, the current density relaxes from the super-critical state 
J > Jc towards the sub-critical state J < J^- The power law behaviour in this regime can be 
attributed to the self-organization of magnetic flux around the critical current density J — J^. 
The self-organization of vortices is driven by inter- vortex repulsion and vortex-pinning center 
attraction alone and is independent of the temperature. This explains the plateau in S{T) 
in Fig. 4, and is also consistent with the apparent universality of S{T) repeatedly observed in 
HTSC §. Moreover, in HTSC the critical current density Jc is 1-2 orders of magnitude less 
than that of low- Tc superconductors over a wide region oi H — T phase diagram. Since the 
relaxation towards Mq oc Jc is a power law, the time required to set up critical state in HTSC is 
much longer compared to low- Tc superconductors. This implies that relaxation measurement 
in HTSC is influenced by the self-organization for much longer time and could be the source 
of non- logarithmic behaviour, ft is important to note that the \ogt relaxation due to thermal 
activation becomes a dominant process only in the sub-critical state. 

The results can be summarized by a t — T plot for the magnetic relaxation in hard type-II 
superconductors which is shown in Fig. 5. Emergence of a new time scale at a temperature 
at which history dependence sets in defines dynamically the irreversibility temperature for 
hard superconductors. The regime T^r < T is analogous to the vortex liquid (VL) phase with 
critical current density Jc = 0. With decreasing T below Tcr, the flux motion becomes rapidly 
viscous as evident from the power law increase in Tq. The overall behaviour is analogous to 
the "supercooled" state in glass formers |l^. In these system, relaxation of the correlation in 
density fluctuation shows scaling behaviour in the /3-relaxation regime jlH which is remarkably 
similar to Fig. 3 here. Relaxation data in HTSC need to be reanalyzed to observe the T 
dependence of the characteristic time scale across the irreversibility temperature and glass 
transition temperature as observed here. Also, in fabricated JJA, SQUID parameter (3j — Xj^ 
of 30 have been achieved |2^ which falls within the parameter range in which the simulation 
results can be applied. 

In conclusion, simulation of the magnetic relaxation in a model of hard superconductor 
shows emergence of a new time scale below the irreversibility temperature. Self-organization of 
the magnetic flux around Jc in this time scale leads to power law decay of the magnetisation. 
Divergence of time scales at a transition temperature leads to a state with finite remanent 
magnetisation (and hence, persistent current). The temperature dependence of the normalized 
relaxation rate is in good agreement with experiments on HTSC. 

Acknowledgements : The author is grateful to D. Dhar and A. K. Grover for useful 
discussions and critical reading of the manuscript. 
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Figure Captions. - 

Fig.l Relaxation of the remanent magnetisation M(r) on logr scale for Aj = 0.1, / = 5 and 
array size 16 x 16. The curves (from left to right) are for temperature T=0.9, 0.6, 0.4, 
0.3, 0.24, and 0.20 to 0.02 in steps of 0.02. Inset: Plot of x(T) for / = 5 in FC and ZFC 
state. The irreversibility temperature is also the crossover temperature Tcr (marked in 
the inset) at which M{t) develops a kink. Also shown is a typical 5x5 array with 
4'r,x, 4'r.y and /r, for a single cell. 

Fig. 2 Log-log plot of T^'^^{T) against T — Tgc for all values of T for which the simulation was 
performed. The symbols V and o are for / = 2, / = 5, respectively. The errors arc less 
than the symbol size. The fitting parameter Tg^ and 7 is also given in the plot. Inset: 
(a) Arrhenius plot Ta{T) vs T~^. (b) The stretching parameter a as a function of T. 

Fig.3 M(t) vs rescaled time T/Ta{T) for 0.05 < T < 0.13 (in steps of 0.01) in the /3-relaxation 
regime. The thick dashed line is a fit to M(t) ~ with b = 0.85. Inset : r/3 as a 
function of T — Tgc for 0.05 < T < 0.14 obtained by fitting the relaxation away from the 
plateau. The full line is a fit to ^ {T - Tgc)"^ with ^ = 1.58± 0.2 and T^c = 0.033. 

Fig. 4 The exponent a oc S (normalized relaxation rate) for T < 0.12. Also shown is the S 
obtained from M(t) around the plateau in Fig.l. Inset : The flux distribution $r,/$o 
across a central row of cells in the array for various r at T = 0.07 for f — 5. 

Fig. 5 The t — T plot showing various regimes obtained from the simulation. The dotted line 
is obtained by interpolating Ti3{T) to Tcr- A and B are the power law regimes with 
exponents a and h, respectively. 
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Fig. 4 
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Fig. 5 
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